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Abstract 

The status of our understanding of relativistic sum rules is re- 
viewed. The recent development of new theoretical methods for the 
evaluation of these sum rules offers hope for further advances in this 
challenging field. These new techniques are described, along with a 
discussion of the source of difficulties inherent in such relativistic cal- 
culations. A connection is pointed out between certain sum rules for 
atomic interactions with charged particles and those for interactions 
with photons. 
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1 Introduction 



In the theory of energy deposition from fast, charged particles and from 
photons, sums of quantum mechanical matrix elements weighted by a power 
of energy can play an important role El E] ■ While the matrix elements 
involved depend explicitly on the excited states of the system under consid- 
eration, these sums can often be shown to depend solely on properties of that 
system's ground state. Such a demonstration involves the use of closure, and 
the resulting expressions are known as "sum rules" . In this article, I review 
work on these sum rules for the case of a relativistic target. 

As an illustration of such weighted sums of matrix elements, calculations 
of the average energy transfer in descriptions of scattering processes lead to 
sums of the form, 

S = Y.EnMn\m\ (1) 

n 

where the sum is over all eigenstates, |n) with energy E^, of the Hamiltonian 
describing the target system, A is an operator representing an interaction 
between incident particle and target, and Eni = En — Ei, with Ei the energy 
of the target's initial state For a dipole interaction, A = z, Eq. 
becomes the well-known Thomas-Reiche-Kuhn (TRK) sum, whereas with 
A = e'''-^/^ it is the Bethe sum, with q the momentum transferred to the 
target, and z and r are position operators. A non-relativistic theory of photon 
absorption involves the interaction A = e ■ pe^^'^, where e and hk are the 
polarization and momentum of the absorbed photon, the magnitude of the 
latter being related to the energy transfer by hk = Eni/c. A relativistic 
theory for photoabsorption has a similar form, but the momentum operator 
p is replaced by the Dirac matrix, a. 

In non-relativistic (NR) theories, evaluation of these sums is often quite 
simple. One need only remember the Schrodinger equation, H\n) = En\n), 
and the closure relation, J2n = ^ with I the identity operator, where 

the sum extends over the complete set of states of the Hamiltonian. For a 
single particle bound in a static potential, we have the Schrodinger Hamilto- 



Then, using A = z, the TRK sum rule is 



3 



= -^^Eni{i\z\n){n\z\i) 



n 

— Y,{^\z\n){n\[Hs,z\\€) + {i\[z , H s\\n) {n\z\i)) 

n 
771 

-{z\[z,[Hs,zW) = l. (3) 



h 



h 

In this paper, I will discuss the effects of relativity on this and other sum 
rules. Unlike their NR counterparts, calculations of relativistic sum rules 
are fraught with difficulties. These difficulties can generally be traced to 
the presence of negative energy states, or from another point of view, the 
existence of positrons. Let us see what happens if we simply replace Hs in 
Eq. Q by the Dirac Hamiltonian, 



Hd = \J3m& + ca-p + Vy (4) 

Then, 



S, = j^^\[z,[Hr>^zm=Q, (5) 

a result of the fact that i^^ is linear in the momentum operator. There is 
obviously no way that this result can reduce to that obtained in the NR 
calculation, indicating that there is something seriously wrong. It is not dif- 
ficult to show that transitions into final states \n) of negative energy have 
precisely canceled transitions into those with energies that are positive |3] . In 
a relativistic theory and assuming the system is initially in its ground state, 
the negative energy states are all filled ("Dirac sea"), so that transitions 
into these states are specifically forbidden by Pauli's exclusion principle. Of 
course, our naive approach utilizing closure in the Dirac representation au- 
tomatically includes these forbidden transitions. To obtain sensible results 
for relativistic sum rules, a method must be found to impose the exclusion 
principle by projecting out the negative energy states. 

The case of photoabsorption has additional difficulties, since the operator 
A appearing in the matrix elements then includes an exponential of Eni- One 
might again try eliminating the explicit dependence on En in favor of if, by 
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but because of the fact that H and z do not commute, this relation is not 
even correct. Therefore, as will be discussed below, calculations of photoab- 
sorption sum rules have resorted to expansions of the exponential in powers 
of EniZ, under the assumption of long-wavelengths, and keeping only the first 
few terms in the expansion. 

The outline of the paper is as follows: In the next section, I discuss how 
studies of interactions between atoms and charged particles or photons lead to 
a consideration of sum rules. In Section I review early efforts to evaluate 
these sum rules using a relativistic model of the target. Following this in 
Section |1[ I present the results of my own work on this subject, studying 
the relativistic Bethe sum rule in the context of stopping power or mean 
energy loss from charged projectiles. I also demonstrate a close relationship 
between other sums that arise in the context of charged particle scattering 
and those related to studies of photon absorption. In Section I consider 
relativistic sum rules from the perspective of Dirac's hole theory, and offer 
some thoughts on how the infinities present in this theory - of energy and of 
the number of electrons - lead to difficulties similar to those described above 
in the context of a single-particle theory. Finally, a conclusion is given in 
Section IHl 

2 Origin of sum rules 

In this section, I review basic ideas concerning interactions between atoms 
and either charged particles or photons, and show how sum rules arise in 
such contexts. I begin in Sec. 12.11 with a discussion of the interaction of 
fast, charged particles with matter, and then in Sec. 12.21 describe the case of 
photoabsorption. 

2.1 Scattering of fast, charged particles 

In the first-Born approximation, the differential scattering cross-section (DCS) 
for an excitation of the target to a state of energy by a fast, charged pro- 
jectile may be written in the Coulomb gauge as ^ 
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[Q (1 + g/2mc2) - ElJ2mc^f 

Here, v is the initial velocity, and z the atomic number, of the projectile, and 
Q = \/m^d^ + q^c^ — m(? is known as the recoil energy. The component of 
(3 = v/c perpendicular to q is denoted by (3t. The quantities, F„ and Gn in 
Eq. ((7j), are defined by 

Fn = E(n|e*'^-^/'*|2), (8) 
i=i 

Gn = E(n|e^^"'-^-^/"«,|z), (9) 
i=i 

where Z is the atomic number of the target atoms. The first term in Eq. (jZj) 
is referred to as the longitudinal interaction and may readily be derived using 
the static, unretarded Coulomb interaction. The other term, that containing 
Gni corresponds to the exchange of a single, transverse photon and is referred 
to as the transverse interaction. 

Use of Eq. (|7j) leads to consideration of sums in various contexts. An 
important example is the calculation of the stopping power, or mean energy 
loss per unit pathlength from the projectile [HllZl- This is found from dan as 

dE — ^ I „ I" , , . 

-^ = T.FniJ dan, (10) 

and the prime on the summation symbol means the sum is restricted by 
energy conservation, such that there is a maximum energy transfer, Eni < 
Em ax- From Eq. (0), it is seen that the above equation involves an integral 
over the recoil energy. It turns out from kinematic considerations that the 
range of this integration is Qm < Q < Quax, where the minimum is related 
to the energy transfer as Qm = E^j2mv^, and Qmux = 2mt>^(M/m)^, with 
M the projectile mass [H Since the sum does not include a complete 
set of states of the target system, closure may not immediately be used to 
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perform the summation. Nonetheless, it has been argued [3] that it is a 
good approximation to extend this sum to include all energies, J2' ~* J2ny 
as long as one also replaces the upper limit of the Q-integral with Qm = 
2mf^/(l — v'^/(?). The errors incurred by this replacement are part of the 
much studied shell corrections [SI 13 CHI HI] • 

When f <^ c, only the longitudinal interaction is important, and then 



dE 

dx 



Eni 



\Fnm 



;i + Q/m^)dQ. 



'Ill 



, Q''{l + Q/2mc^f 

It is advantageous to perform the summation before the integration, but one 
must first account for the fact that the lower limit of integration depends 
explicitly on En- A common approach is to split the Q-integral into two parts: 
the low- and high-Q regimes, separated by a cutoff, Qc- Then the integration 
limits for the high-Q regime do not depend on the energy of the excited state, 
and the order of summation and integration may be interchanged in the 
expressions for these integrated quantities. Thus, we are led to consideration 
of the Bethe sum jT2j , 



iq-fj/h 







;i2) 



Si{Q) = Y.Eni\Fnm^ = Y.Y.Em me 

n j=l n 

where I have followed Bethe in neglecting the small correlation terms involv- 
ing two different electrons. 

For the low-Q region, the integration must be done first. To accomplish 
this, a long-wavelength approximation is introduced by expanding the expo- 
nential in Fn{q) as e*'^'^-' ~ 1 + iq-fj. Then, neglecting terms 0{Q /mc?) gives 
-^n(^) ~ Qfn/Eni, with the optical dipole oscillator strength defined as 



fn 



'2mEni 



Y,{n\zj\i) 



(13) 



(14) 



This leads to the TRK sum rule, which for a system of Z electrons is 

S^=J2U = lim—Si{Q) = Z. 

For a relativistic projectile, the transverse part of the DCS must also be 
included. The appearance of Eni in the denominator of this term complicates 
its treatment in terms of sum rules. I have recently shown how to treat the 
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high-Q case by expanding the denominator in powers of the small quantity 
{Q — Eni)/Q [ini, leading to sums of the form 



z 




(15) 



j=l n 



For the low-Q case, expansion of the denominator is not appropriate, 
since Q — Eni and Q may be comparable in magnitude. When the target 
system is non-relativistic, Fano showed that the transverse contribution is 
proportional to 5*1 for small Q. In the case of a relativistic target, I have also 
shown [13j that Fano's approach must be generalized, as will be discussed in 
Section lOl 

2.2 Photon interactions 

The study of sum rules for photon-matter interactions grew out of a conjec- 
ture by Gell-Mann, Goldberger, and Thirring (GGT) ^^l- In the course of 
investigating the implications of causality on quantum theory and, in partic- 
ular, on scattering amplitudes derived from field theory, these authors arrived 
at a dispersion relation between the dispersive and absorptive parts of the 
photon forward scattering amplitude. Let us now look at how this leads to 
sums that are quite similar to those discussed in the previous section. 

In their paper, GGT considered the matrix element of the 5* matrix, 
S = 1— iR, for scattering of a photon of four- momentum k" and polarization 
a into one of four- momentum k' and polarization A, by a matter system 
which is ultimately left with its energy unchanged. They wrote. 



{Xk'f\R\iak") = -I dV / d''x"{Xk'f\V{x')V{x")\iak")T^{t' -t") 



with the step function ri{t' — t") vanishing unless t' > t", in which case it is 
unity. The interaction potential is 




(16) 




(17) 



where the current operator j and the transverse vector potential. 
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•' s=l ^ 

are operators in the interaction representation. The operator creates 

(annihilates) a photon with momentum hk, frequency u = kc, and polariza- 
tion es perpendicular to k. 

Since the initial and final states each have a single photon present and V 
creates or annihilates one photon, it is clear that in Eq. ()16|). there are two 
possible ways for the system to pass from the initial to the final state. Either 
a photon is first emitted and then one is later absorbed, or vice-versa. This 
yields two separate terms in Eq. from which the scattering amplitude 
becomes |JJJ 



{f\jx{x')\n){n\j^{x"M 



^~ik-x' +ik' -x" 



^ {fljajx') \n) {n\jx{x") \i) ^k-x'-ik'-x" 
Ei — huj — En + ie 



Ei + hu — En + ie 

(19) 



By considering the analytic properties of this quantity, GGT derived the 
following dispersion relation between the real (TZe) or dispersive and the 
imaginary (Xm) or absorptive parts of the forward scattering amplitude for 
photons of frequency u: 



TZefiu) - 7^e/(0) = P / du 

TT Jo 



= ^Pfd.'^. (20) 

where f{uj) = F^)^{k,k), P denotes the principal value of the integral, and 
use has been made of the optical theorem, Zmfioj) = Luaabs{i^)/4:TT, with 
c"a6s(<^) the total absorption cross section for a photon of frequency u. For 
scattering from bound electrons, TZef{0) = (Rayleigh scattering). Then, 
with the conjecture that for u oo, f{uj) approaches its free electron value, 

7^e/(oo) = -e^/m, (21) 

GGT obtain the sum rule, 
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/ ducrabs{uj) = . (22) 

JO m 

As I will now demonstrate, this is closely related to the TRK sum rule, 

Eq. Q. For a single non-relativistic electron, we have 

m = ^{p,S{x-f)} (23) 

Then the total photoabsorption cross section may be found directly from 
Eq. ()19|) and the optical theorem, and is given by 

CTaU^) = Y.\(^\P-^ - Em) + + i^n^)], (24) 

where I have taken the polarization A along the x direction and the momen- 
tum k along z. If the system is initially in its ground state, Eni > 0, and 
the second delta function will not contribute for positive frequencies. Then, 
expanding the exponential, we obtain 



f°° J / N 4:7i'^e'^h^ 1 /,. , ....2 Eli |.\i2 
/ duaabsiuj) = —-^2^— [\{n\p^\t)\ + \{n\p^z\t)\ 

10 fn ^ -Ly-ni \ He 



-^^Mp^\i){n\p^z'\i)+H.c.]^, (25) 

with H.c. the hermitian conjugate. The first term in parentheses is the dipole 
term El., the next is the quadrupole E2, and the remaining terms represent 
retardation effects on the dipole interaction, denoted as EV^^. 

In the dipole approximation only the first term is kept, and using the 
identity Px = im[H,x]/h, we see that the GGT sum rule, Eq. (j221), becomes 

^j:Em\{n\x\z)\' = l, (26) 

which is just the dipole oscillator strength sum rule of TRK, as claimed. The 
important point is that Eq. (j^^ is much more general in that it has been 
derived without the assumptions of long wavelengths, a dipole interaction, 
and a non-relativistic matter system. Unfortunately, as we will discuss in the 
following section, the conjecture Eq. ()21|) is incorrect and so therefore is the 
GGT sum rule. 
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3 Review of early work 



The conjecture of GGT was based on the notion that for very high energy 
photons, one might reasonably suppose that the binding energy of the elec- 
tron will have little import, and can therefore be neglected. Nonetheless, 
GGT expressed concern about the correctness of this conjecture. In partic- 
ular, they pointed out that the photoelectric cross section was then believed 
to have a 1/a; dependence at high frequencies, leading to a logarithmic oj- 
dependence on the right-hand side of Eq. (j^ . 

The apparent divergence of TZef{uj) at large uj can be circumvented by in- 
troducing a procedure that was originally suggested by Brown jTH] • This pro- 
cedure involves defining the relevant cross section as the difference between 
the photoelectric cross section and (Tpp^, the cross section for pair production 
in which the electron of the pair is produced in the occupied Is state. As was 
later discussed in detail by Erber ^^1; this definition of the cross section cor- 
responds to scattering by the bound electron, alone. That is, it is in fact the 
difference between the cross section for a one-electron atom and that for an 
isolated proton Erber than argues that the logarithmic dependence of 
the dispersive amplitude that was mentioned above is exactly canceled by this 
subtraction procedure. It should be noted that along with this subtraction of 
cross sections, one must include a corresponding subtraction of forward scat- 
tering amplitudes, which has also been carefully discussed by Erber |TS|. As 
is shown in his Eq. (3.1), the correct amplitude is not then obtained directly 
from the scattering amplitude of GGT, our Eq. (fT^ . Rather, by carefully 
considering the various processes involved in scattering from the atom and 
also from the free proton, and in particular, the processes involving creation 
and annihilation of electron-positron pairs, Erber has argued that the correct 
amplitude requires the replacement of ie by —ie in the second term of our 
Eq. (fT^ (compare Erber 's Eq. (3.1)). 

Shortly after publication of the paper by GGT, Levinger and co-workers, 
using the Brown/Erber procedure, showed that the conjecture, Eq. (pT|). 
still turns out to be incorrect jUHTHIHl. These authors began [T7[[TH] with 
numerical studies of the dipole oscillator strengths for a single electron in the 
field of a Ph nucleus undergoing transitions from the Is into higher discrete 
states. They found that the retarded, relativistic oscillator strengths for these 
discrete transitions were only about 80% of their non- relativistic values ^7] . 
For transitions to the continuum states, they obtained values from the work 
of Brown OU] and Hulme [21]. Including Jacobsohn's non-relativistic 



11 



calculations of the electric quadrupole transitions {E2) '1^, they found a 
summed oscillator strength of approximately 0.86, in stark contrast to the 
value of unity predicted by GGT (inclusion of retardation in E2 changes this 
value to 0.87 |4|). They argued that this result should be accurate to within 
a few percent, even considering the neglect of higher multipoles. 

The following year, Levinger, Rustgi, and Okamoto addressed this ques- 
tion through the evaluation of sum rules To overcome the difficulties 
arising from the existence of negative energy states, these authors used 
a semi-relativistic Hamiltonian, which does not have the negative energy 
eigenstates of the original Dirac Hamiltonian. The important difference be- 
tween their Hamiltonian and Hs is the inclusion of the mass-variation term, 
Hm = —p^/8m^c^ (other correction terms in their Hamiltonian are at most 
linear in the momentum operator, and therefore do not contribute to the 
sum rule). Since the negative energy states have already been removed, the 
Brown/Erber subtraction procedure discussed above is not appropriate here 
(nor is it needed in the related sum rule calculations discussed below), since 
that procedure is specifically related to pair production, which is not included 
in this approach. 

For the unretarded El term, Hm adds a relativistic correction of (see 
Eq. ©) 



El 



rel 



1 



(27) 



for an isotropic initial state, with To = {i\p'^ / 2m\i) . Using this same Hamil- 
tonian, the E2 contribution was calculated relativistically, as were the retar- 
dation effects, EV^^, on the dipole transitions. They found 



and 



E2 



1 



Ibmc^ 



+ 4r— 



dr"^ 



dr 



(28) 



El 



ret 



1 



15mc^ 



' nd^V , dV" 
r —-TT + 4r— 

ctr^ dr , 



(29) 



The terms involving V in E2 are exactly canceled by EV^*. This is one 
illustration of the well-known fact that higher multipoles and retardation 
each contribute at the same order of magnitude as relativistic effects, and 
that it is inconsistent to include one without the others 
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Combining results, these authors obtained 



/(oo) = --(l-^)^ (30) 

for the GGT sum rule. They pointed out that the presence of Tq in the 
denominator represents the increase of mass due to the electron's kinetic en- 
ergy, consistent with the fact that its presence results directly from inclusion 
of Hm in the Hamiltonian. To obtain a numerical result, they used the value, 
To/mc^ = = 0.18 for Pb, with a the fine structure constant. This 

yields a summed oscillator strength of 0.82, in satisfactory agreement with 
their earlier result of 0.87 JH]- They observed that exact agreement should 
not be expected, due to the approximate nature of their numerical calcula- 
tions f7|^], as well as the neglect of contributions OiZ'^a^^ in the sum rule 
calculations. 

The above result was confirmed a decade later by Goldberger and Low 
|25j . who calculated the forward scattering amplitude in the limit of very 
high photon energies. They obtained the result, 

/(oo)^--(l-^ + ^ + --), (31) 

which is in agreement with Eq. (j30|l to lowest order. An exact expression for 
/(oo) in terms of quadratures was also derived j2Sl- 

Over the next several years, a number of papers considered the integrated 
photoabsorption cross section [2H1 123 I2HI 1211 EDI , using analytical methods 
based on closure and generalizing the sum rules by: (a) including a world- 
scalar potential and utilizing the Foldy-Wouthuysen transformation (FWT) 
[HT] to obtain a semi-relativistic Hamiltonian ^7] , and (b) studying a many- 
electron system with interactions . In spite of Levinger's work jlj El QB] ) 
Gerasimov concluded that the GGT sum rule, Eq. (j221), was in fact correct 
[21]. By studying specific models, Matsuura [28j obtained results in agree- 
ment with Levinger's 0], and pointed out that Gerasimov had incorrectly 
assumed that j[iS) is analytic everywhere except on the real lo axis, invali- 
dating his conclusions. 

In 1985, Schmitt and Arenhovel [32^ studied the GGT sum, using both 
analytical and numerical techniques. For the analytical efforts, they studied 
the leading relativistic corrections to the sum rule in three different models 
of the target particle: (a) non-relativistic (NR) or Schrodinger, (b) Foldy- 
Wouthuysen (FW) using the FWT j2T| to obtain a semi-relativistic model. 
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and (c) Dirac using Hj:, and introducing projection operators to exclude 
contributions from the negative energy states. In the latter case, they used 
an expansion of the projection operators in powers of Hu/mc^, so those 
results are also semi-relativistic. Their use of the FWT to transform the 
multipole operators for the FW model represented a technical improvement 
over previous efforts |31 |2Z1 , wherein only the Hamiltonian and states had 
been transformed. Nonetheless, all their results were in agreement with that 
of the earlier work. 

These authors also numerically calculated f{oo), using the exact (fully 
relativistic) expression of Goldberger and Low (see Eq. (4.1) of |2H]) for a 
particle in a square well. Their results were such that /(oo) and the FW 
calculation of the GGT sum differed by an insignificant amount (less than 
0.1%). The implications of this result for our understanding of relativistic ef- 
fects is not completely clear, however, since their chosen binding energies are 
rather small, being at most about 5% of the particle's rest energy. Nonethe- 
less, they did find a small, but significant, contribution to the GGT sum from 
the pair production term, apps- 

Shortly thereafter, Leung, Rustgi, and Long [SB] considered relativistic 
corrections to the Bethe sum rule. In a non-relativistic theory, this sum can 
be readily evaluated exactly including all multipoles: 



(32) 



where I have again assumed an isotropic initial state. To obtain the rela- 
tivistic corrections, Leung, et. al. [23] used an approach similar to that of 
|27|, performing the FWT on Hd for a one-electron system, obtaining 



^iW)=Q(i-^(^Ip'I^))> (33) 

where \i) is the transformed initial state, and terms 0{Q^) and 0{Q'^p^) were 
neglected. Leung later used this result to obtain the leading relativistic 
corrections to the Bethe stopping power. Note that this result is obtained by 
transforming the Hamiltonian and states, alone, while leaving the operator 
Q^q-r/h unchanged. Thus, this procedure is inconsistent in its treatment of the 
various operators. We can also see that the result is incorrect by comparing 
to that for a relativistic free-particle at rest [H], 
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The result of Eq. (jH^ should reduce to this expression if the expectation value 
of is set to zero, and is thus seen to be missing the term proportional to 
Q^. A proper treatment using the FWT and transforming all operators was 
later formulated [HS], obtaining 

S,{Q) = g (l - -^2^\P^\^) - QV'^rnc^) , (35) 

in agreement with Eq. (j34p . Other sums, involving various moments of the 
energy transfer, were also evaluated 

It was suggested in that use of a no-pair Hamiltonian and/or 
field-theoretic techniques might be useful for consideration of many-electron 
atoms. Such techniques were later used by Aucar, Oddershede, and Sabin 
jSTj in their work on the relativistic corrections to the Bethe sum. They in- 
troduced a second-quantized no-pair Hamiltonian and projected the operator 
^iq-Tj/h Qj^^Q space of positive energy states, thereby excluding transitions 
into the negative energy states. These authors mention that it is straightfor- 
ward to use their approach for the case of a many-electron atom, but they 
explicitly consider only the one-electron case. Unfortunately, although the 
techniques used in this approach are formidable, it does not appear to the 
present author that there is any simple way to obtain detailed results from 
their expressions, which are still in the form of a sum over a complete set of 
eigenstates of the Hamiltonian. 

4 Recent advances 

4.1 Relativistic Bethe sum rule 

In this section, I describe the work I have recently done on the relativistic 
Bethe sum rule. I will discuss the FWT jSI] , derive a formally exact expres- 
sion for the sum rule, and then explain why the results of [321 ESI only 
valid when q is small. I will also show how to obtain an exact treatment of 
the g-dependence of the Bethe sum rule, using an FWT-type transformation 
introduced by Douglas and KroU ||38|- When combined with the evaluation 
of sum rules for the transverse part of the DCS, these results permit the 



15 



calculation of relativistic corrections to the stopping power that are valid 
for a fully relativistic projectile [12]. At the end of this section, I outline a 
method that appears to give highly accurate expressions for the TRK and 
Bethe sums even for very strongly bound systems. 

To implement the FWT, one inserts the identity operator in the form 
of e~^^e^^ , with U chosen to be Hermitian so that e*^ is unitary. For this 
preliminary discussion, I will assume that U is the exact FWT operator for 
a general many-electron Hamiltonian H. Then, for example, 

{n\He''^-''^/^\i) = {n\e-'^e'^He-'^e'^e''^-'^^/^e-'^e'^\i). (36) 

Thus, e*^ induces a canonical transformation of the states and operators, 
with the identification, \n) = e*^|n) and H = e*^ife~*^. Then, defining 
(f-dependent operators through the replacement of pj by pj + q due to a 
transformation of the form, ^A ^ Ai{q) = e~^^'^^^^A4e^^'^^^^, we have 

{n\He"^-'''/^\i) = {n\e'^-''^'^H{q) (f^^^'^e''^) |z). (37) 

The FWT operator, e*^, is chosen so as to bring the Hamiltonian into a 
"block diagonal" form, explicitly separating the positive and negative energy 
subspaces. Now these positive energy states are still 4-component spinors, 
but their third and fourth components vanish. Therefore, though e^^^^e'"^^ 
is a 4 X 4 operator, only its upper-left 2x2 block will contribute to the matrix 
elements in our sum. Defining this block to be Rj, we then have, 

{n\[H, \i) = {n\e"l''^'^{H{q)Rj - RjH) \i), (38) 

and 

{i\e-''^-''^/^\n) = {i\R]e-''^-''^/^\n), (39) 

which when combined and summed, yield a formally exact expression JH] 
for the Bethe sum, Eq. (O, 

SiiQ) = j:m] {H{q)R, - R,h) I?). (40) 
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To this point, I have assumed that U is the exact FWT operator. How- 
ever, given that this operator is not known for the case of bound electrons, it 
will eventually be necessary to introduce approximations for e*^. Nonethe- 
less, the fact that Eq. ()40|) is an exact result for the sum rule implies that 
the excited states have been treated exactly, and an approximate treatment 
of the operators appearing in the resultant expectation value does not alter 
this conclusion. In particular, this expression remains valid even when a 
semi-relativistic approximation to the FWT is used under circumstances for 
which the original summation includes significant contributions from highly 
relativistic states |n), such as with very large values of Q. 

Let us now see why the approximations used in [331 EH] for a single electron 
yield results that are only valid when q is small. These authors used the 
approximation of |31j . 

gic/ _ ...gic/3gic/2gic/i^ (41) 

which involves an expansion in powers of the small quantities p/mc and 
V/mc^, with ~ 0{l/m^). The calculation of the FW Hamiltonian jBT] 
then proceeds by an expansion in powers of the 11^. For example, the first 
step of the FWT gives 

H, = e'^^i/e-*^^ = H + i[Uu H] - ^[U,,[U^, H]] + ■ ■ ■ , (42) 

which is an expansion in powers of p/mc, since Ui = —i(3d ■ pjlmc. Higher 
order corrections to the Hamiltonian are obtained from and Vj, in a simi- 
lar fashion ^31j. Now, Hi{q) is obtained from Hi by the replacement of p by 
p + q. This means that the above expansion in powers of p/mc leads directly 
to an expansion in powers of q/mc. Since this approach has started out with 
a p-expansion, it ends up with a g-expansion and is therefore incapable of 
handling large values of q. Therefore, unlike the NR Bethe sum rule, the 
relativistic corrections of jSSl do not include contributions from all mul- 
tipoles. In fact, the difference between the free and bound particle results, 
given in Eqs. (j34j) and (j35p . arises only from the dipole contributions. There- 
fore, the latter equation gives us little more information than the corrections 
to the TRK sum rule. 

An expansion in powers of the potential energy V does not lead to the 
difficulty described in the previous paragraph, at least as long as V commutes 
with e'''-''/^ so that e-^^-^^/^^e^^"-^/^ = V. This indicates that the way to 



17 



circumvent the restriction to low q is to avoid, or at least delay, the expansion 
in powers of p. Douglas and KroU [23) have introduced a modified FWT that 
only expands in powers of V [3H1 EH HH] and is therefore suitable for this 
purpose. This transformation (DKT), which is still of the form of Eq. ()41|) . 
takes Ui = Uf for a many-electron system, where 

e'^f = f[e'^\ (43) 

k=l 

and the free-particle FW operator for the k^'^ electron is given by 



^±iu, ^ 1 ^ ^^2 ^p^^^.p^^^ (44) 

'2efc(efc + mc2) 



with efc = yplc^ + rn^d^. 

For if, I have used a no-pair many-electron Hamiltonian [36j 




ii = Vi/^ + /:+ Vi4J/:+, (45) 



where is the single-electron Dirac Hamiltonian for electron fc, Vfc; is the 
interaction (Coulomb-Breit, for example) between electrons k and Z, and 
^+ = nfc=i'^+(^) with £+(/c) the projection operator onto the space of 
positive energy states of It is true that the Vki part of the potential 
energy does not commute with e^^'^^^^^ so it is not precisely correct to claim 
that a l^-expansion will yield the exact g-dependence for this Hamiltonian. 
Nevertheless, the contributions to the sum rule from Vki (apart from its 
effect on the eigenstate which will not depend on q) turn out to be quite 
small fSJE], so the error incurred by omitting this g-dependence should be 
neghgible. 

The first step of the DKT transforms H into 

_ z 

H = Y.eu + V. (46) 

where I have written V = e'^f(J2f=iVk + C+J2i>iVki'C+)e-'^f to represent 
the potential energy terms. As I will discuss below, these terms appear to 
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lead to very small corrections, so will not be further considered in any explicit 
manner. Using this as our transformed Hamiltonian, we have 



H{q)=jZet + eM + V. (47) 

Furthermore, 

giC//(9) = gic/,(5l -Qgit/.^ (48) 

and 



^UM ^ 1 =(6,(g) + mc^ ± ■ (p, + q)). (49) 

All these expressions may now be incorporated into Eq. pn|) and ex- 
panded in powers of pj, while keeping the exact dependence on q. When this 
is done, one finds the relativistic Bethe sum rule to be j3T] 

Sr{Q) = ZS',{Q){1 - AjiQ) - A,,(Q)), (50) 

with 



A/(Q) 



Sm^c^Z (1 + Q/mc^y 



Sm^c^Z (1 + g/mc2 



E(obllo), 



(51) 



and 



^'•(« - - 2,n3c^Z(iyQ)f/^)^(lTQ/2,„c^) g<°l^^ ' ^''^ ^ 

Note that in A^,(Q), the coefficient of V^Vj- vanishes identically. This is 
somewhat surprising, given the complicated nature of the calculations leading 
to this result. Expectation values of the remaining spin-orbit operator are 
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usually smaller then some of the two-body terms that arise from Vm, and 
which I have omitted jini SI] because they are known to be quite small 
|12l 113 El] ■ Therefore, to a high level of accuracy we have 

A,(g) = 0, (53) 

and the potential energy of the electrons does not appear explicitly in the re- 
sult, entering only through its effects on the ground-state wavefunction of the 
target system. These results give the relativistic corrections to 0{p'^/m'^c^) 
and (9(p^K/m^c^). Since they are valid to all orders in Q/mc^, they include 
contributions from all multipoles, unlike earlier studies EH]- It should 
also be remembered that the Bethe sum arises from the longitudinal part 
of the DCS and as such, does not involve retardation, in contrast to the 
photoabsorption sum of GGT. 

The TRK sum rule may be obtained directly from the above results as 
Si = Z(l- A/), with im 



% = Jj^pipm - ^hom (54) 

The first term is just that found by numerous others [U 123 EHl IIH ES], but 

the second term had not previously been known. The reader may recall 
that the TRK sum rule for a single electron in the field of a Pb nucleus 
was evaluated (along with electric quadrupole and retardation contributions) 
using perturbative methods in |1]. It is consistent with these methods to 
evaluate the term using an NR wavefunction. For the one-electron case 
and a Pb nucleus, this gives a contribution of — 35(Za)^/8 ~ —0.53, whereas 
the first-order term is 0.30. Thus, it is quite clear that this approach fails 
badly for such strongly bound systems. For neutral atoms, on the other hand, 
we may expect that the average over electrons will diminish the errors and 
allow the perturbative approach to give a reasonably close approximation. 
This notion is supported by the results 1 obtained in [H]. Nonetheless, it 
needs to be checked with the more accurate approach described below, using 
DKT wavef unctions. 

Numerical studies of Eqs. and indicate that perturbative calcu- 
lations of these sum rules for neutral atomic systems yield accurate results 
over a wide range of Z |41j. Further work is needed, however, when Z > 70 
and Q is small, if one wishes to achieve an accuracy of 0.5% or better ^T] . 
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To obtain this greater accuracy for high Z, it is possible to use the methods 
already described, but to forego completely the expansion in powers of pj. 
That is, we may insert Eqs. (gH), (gUI), (071), and (gHl), directly into Eq. (gUl), 
and evaluate the expectation value using the accurate DKT ground state of 
H. Expectation values using these states, and carefully including the effects 
of Vki, have been used extensively by Hess and co-workers [SnillSllini- The 
success these authors have had with such calculations, keeping only low-order 
terms in V, makes it quite likely that this approach will also yield highly ac- 
curate results for the sum rules. 

For the 0{V) terms calculated without the p-expansion, numerical eval- 
uation of the Bethe sum for nonzero Q appears to offer quite a significant 
challenge jl7j. Therefore, I have begun with consideration of the TRK sum, 
5*1 = Z{1 — Af — Ay), keeping only 0(y) in Ay. Now, A/ has the relatively 
simple form. 



here. Though accurate numerical calculations of Aj and Ay are yet to be 
completed, I have obtained order-of-magnitude estimates of these quantities 
for a single-electron system by using NR hydrogenic wavef unctions. The 
results are intriguing, and they will be the basis of the discussion in the 
remainder of this section, which will mainly be speculative in nature pending 
the outcome of work in progress. 

For all values of Z, these rough estimates indicate that Ay is much smaller 
than Af. Considering once again a single electron in the field of a Pb nucleus, 
I find Af = 0.17, while Ay = 0.002 . If instead we have a nucleus with 
Z = 120, then Af = 0.27 and Ay = 0.008. We see that even for such extreme 
conditions, Ay/Aj a; 0.03, and Ay itself is only about 1% of the total sum. 
For most physical systems, only a small fraction of the electrons will have 
binding energies approaching those just considered. For neutral atoms, I 
have found A/ to be no more than about 0.03 in the perturbative approach 

1^ . If the relative size of the two terms remains even approximately the 
same, than we will be justified in neglecting Ay. If this is confirmed by the 
more accurate numerical methods using many-electron DKT wavef unctions, 
then Eq. ()55|) will give us the effectively exact relativistic corrections, 
to the TRK sum rule. Since this expression does not involve the potential 




(55) 



but the expression for Ay is quite complicated, so will not be written out 
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energy V and is rather simple in form, it should be relatively easy to evaluate 
even with realistic wavef unctions. 

We may hope for more. In studies of the Bethe sum, it was found that 
the relativistic corrections were never much larger than their value at Q = 
|41j . Therefore, given the relationship between the TRK sum and the Q —>■ 
limit of the Bethe sum (see Eq. (HD)), it seems probable that Av{Q) will be 
negligible over the whole range of Q. If this can be demonstrated, then the 
corrections to the Bethe sum rule will be given to high accuracy by setting 
V = in the operators appearing in the expectation value of Eq. fl40|) . This 
reduces to the simple expression [T31 ESI , 



Si{Q) 



2 Z 



E 



Q +q-Pj ^ q-Pj 







(56) 



Then this approach, using the DKT without expanding in powers of p, could 
be used to obtain accurate stopping powers for targets consisting of even the 
heaviest elements. 



4.2 Photoabsorption sum rules for scattering of charged 
particles 

In this section, I will discuss the summations that arise from considerations 
of the transverse part of the DCS, given in Eq. ((Zj), and show they are closely 
related to the GGT sum. In [13] , I have evaluated sum rules for these contri- 
butions to the stopping power using both low- and high-Q approximations. 
In the high-Q region, the recoil energy and energy transfer cannot differ sig- 
nificantly PlEl, so it is a good approximation to make an expansion in powers 
of {Q — Eni)/Q, which leads to sums of the form given in Eq. (fT3j). Although 
these sums involve the Dirac operators aj, they can be handled with the same 
techniques described in the previous section for treatment of the longitudinal 
contributions. Details of this approach, along with corresponding results for 
the stopping power, are given in [13,. 

The low-Q transverse contributions to the stopping power are given by 



|A-G„| 



[Q (1 + g/2mc2) - Elj2mc^Y 



(l + Q/mc^) dQ. 



(57) 
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Fano's treatment of Bt began with the change of variables [S], 

QmiQm + 2mc^) 



cos^ ip 



Q{Q + 2mc2) 



(58) 



giving q = Eni/v cos ip and 



ZE.^ 



sin^ Tp 



ni cos 



s2t/.e [1 - /^^ COS^ l/j^ 



n e 



iqZj /h 



a 



JX 



^(cos^V^). (59) 



Here, is the angle between P and q so that (3t = (3 sin ^z^, and is the 
value of %l) when Q = Qc- Fano then set g = in the exponential and 
related the matrix element of ax to that of x, reducing this expression to one 
proportional to the dipole oscillator strength, [Sj. 

In obtaining the relativistic corrections to the stopping power, I found it 
necessary to keep terms 0{q^) from the exponential in Eq. (|59p in order that 
the results for the stopping power were independent of the choice of cutoff, 
Qc 53 • Rewriting q in terms of Eni and expanding the exponential allows 
the -^-integral to be done, and yields an expression for Bt involving the two 
sums. 



E 



1 



E ■ 



(60) 



and 



7 = 1 




Z 

E {iWkx\n)(n\zjajx\i) + H.c. 
i,fc=i 



(61) 



where the terms linear in Zj vanish due to conservation of parity. Since 



[H, Xj 



-ihcajx (apart from small two-body terms), Eq. (jHIH) is just the 



relativistic TRK sum, representing the unretarded dipole contributions El. 
Comparison to Eq. (^3)) reveals that Eq. is just the relativistic form of 
the combination, E2 + EV'^''^ , considered in the work on photoabsorption sum 
rules. It is thus seen to have been no accident that for a relativistic treatment 
of the target, I needed to keep higher-order terms in the expansion of the 
exponential of Eq. ()59j) . As has been pointed out previously, the effects of 
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higher multipoles, retardation, and relativity all lead to contributions of the 
same order of magnitude. Keeping only the leading term in the expansion 
excludes the E2 and EV^^ contributions, and cannot yield a consistent rela- 
tivistic formulation. We also see here a direct relationship between the GGT 
sum, which arose out of studies of photoabsorption, and the sums related to 
the transverse contributions to the stopping power. 

Though we have arrived at this conclusion from consideration of the low- 
est order multipoles, it can be shown that this relationship holds for all 
multipoles. If the lower limit of the integral in Eq. is avoided, then we 
may directly extract from Eq. ()59|) the sum. 



which is just the relativistic generalization {p a) of the integrated pho- 
toabsorption cross-section in Eq. Thus, the sums for the transverse 
contributions to the stopping power have contributions that are exactly the 
same as those appearing in the photoabsorption dispersion relation of GGT. 

We may understand the appearance of the photoabsorption sum in the 
context of charged particle scattering in the following way: First, note that 
this contribution to the stopping power represents the exchange of transverse 
photons between projectile and target. If the target is initially in its ground 
state, it cannot emit, and can only absorb. Then, observe that in Eq. ()59|). 
the projectile states have already been integrated out, leaving a description in 
terms of the target states alone. Thus, this equation describes the absorption 
of photons by the target, and it is completely reasonable that it contains the 
photoabsorption sum. 

Nonetheless, the reader may note two important differences between Eqs. 
(j24|) and (j59p . First, for the charged particle case, the energy-momentum 
relation given after Eq. (j58|) is not that usually seen for photons. This is be- 
cause these are virtual photons, for which energy and momentum can take on 
any value, and each may do so independently of the other. One may imagine 
that both projectile and target are constantly emitting and reabsorbing vir- 
tual photons, and occasionally one of these photons gets exchanged between 
the two. The photons that are exchanged must be ones that lead to con- 
servation of both energy and momentum in the overall interaction between 
projectile and target, since that is a real process. This requirement leads 
to the relation, q = Eni/v cos (except at extremely high incident energies. 




2 



(62) 
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approaching the TeV range for protons) instead of Eni/c. 

The second difference has to do with the presence of the T/'-integration in 
Eq. (jSni)- This integral simply reflects the fact that contributions from all 
possible final states of the projectile have already been summed. Therefore, 
unlike in the photoabsorption case where there is only one possible momen- 
tum transfer for each Eni-, here there is a range of momentum transfer, or 
equivalently if), that is possible. One may also note that the upper limit, 
cos^'j/'c, of the 'j/'-integral depends on Eni- Therefore, there are additional 
terms that differ in form from those resembling the photoabsorption sum. 



5 The trouble with relativity 

We have seen that significant challenges arise in the evaluation of sum rules 
for relativistic systems. Unlike their non-relativistic counterparts, the use 
of closure for relativistic sum rules requires the introduction of special tech- 
niques, such as the FWT or projection operators exclude contri- 
butions from negative energy states. Direct use of closure without such 
techniques in a single-particle theory of a hydrogenic system includes con- 
tributions from transitions of the electron into negative energy states. As 
previously stated, these contributions precisely cancel those from transitions 
into positive energy states |3], yielding unphysical vanishing results for the 
sum rules. 

If one views these sums from the many-electron perspective of Dirac's 
hole theory, on the other hand, the negative energy states are filled (Dirac 
sea) and transitions into these states are automatically excluded. Then one 
might expect an absence of such difficulties when using the closure relation, 
= 2r, with \n) now a many-electron state including the filled Dirac 
sea. However, one still obtains a vanishing result when an average over the 
direction of q is taken (the mathematics is essentially the same as for single- 
electron theory). The explanation of this result is no longer straightforward, 
but as I will now describe, it is related to the appearance of infinities from 
the participation of electrons in the Dirac sea. For the following discussion 
let us consider the Bethe sum, 

Si{Q) = Y.llEm\{n\e"i-''^"^\T)'\ , (63) 

j n 

and now includes all negative energy electrons, as well as those with 
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positive energy. For simplicity, I will restrict the discussion to a single, free 
electron at rest in the presence of the Dirac sea. Then with \i) and \n) given 
by suitably symmetrized sets of momentum states, the matrix elements for a 
given q are proportional to 6{pj ' —pj — q), and only one final state contributes 
to the sum over states for each j. For bound electrons, this probability 
distribution will be broadened, as has been discussed in reference to the 
Bethe surface Pfl El EH] , but the basic arguments presented below still hold. 
Therefore, the conclusions should be quite general. 

In consideration of this system, we may identify two possible sources of 
infinity. The first is the infinite energy, Ed, of the Dirac sea contained in 
each of and Ei. The other is the fact that J2j runs over an infinite number 
of electrons. Let us consider each of these in turn. 

For a given j, we have a sum over states \n) and a subtraction of the 
energies, Eni = En — Ei. If the subtraction is taken first, explicitly calcu- 
lating each individual term prior to the summation over states, there is a 
precise cancelation of the infinite E^ between En and Ei. Then a finite re- 
sult is obtained, which is proportional to for small q (this calculation is 
straightforward for the free-particle states presently under consideration). 

On the other hand, when the sum over states is done first using closure, 
we obtain a very different result, 

S^iQ) = Y.{Ae-''-''^^''[H,e''-^^/^]\t) = c^(^|a, ■ gj^). (64) 
j j 

For each j, the result is again finite, but here it is linear in q and therefore, 
upon averaging over the direction of q, we obtain the vanishing result referred 
to above. What has happened is that we have reformulated the sum into a 
difference of two terms, 

n n 

each of which is infinite. While the infinities still cancel, the resulting finite 
result is of a very different form as that found by explicit calculation of 
each individual term, followed by the summation. Because of the presence 
of infinities, it is incorrect to interchange the order of these mathematical 
operations, the subtraction and summation being considered here. 
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One may ask, which method is the correct one? We will argue that it 
is the first one, where the subtraction is taken before the summation, even 
though, as we shall see, this yields an infinite result. 

To begin with, I remind the reader that the Bethe sum arose in a cal- 
culation of the mean energy loss from fast charged particles to matter. The 
definition of this sum thus involves terms, each of which has a factor of the 
energy transfer, Each term to be summed is proportional to this energy 
difference; that is, according to the original derivation of the sum, the differ- 
ence should be taken before the summation over states. As further evidence, 
note that our result should reduce to the NR result in the appropriate limit. 
The NR result is also proportional to and does not vanish upon averaging 
over the direction of q. Therefore, we conclude that the correct result is that 
obtained when the summation over states comes after the subtraction of en- 
ergies. The alternative, where closure is used directly to obtain a sum rule, 
is seen to be faulty, as a result of the interchange of mathematical operations 
in the presence of infinities. 

Next let us consider the sum over electrons. As I have said, for any 
given free electron only a single transition contributes to the sum, that for 
which ' = Pj + Q- For a positive energy electron initially at rest, the en- 
ergy difference between these states is just Q < Qm- Transitions involving 
energies greater than Qm cannot occur. In this sense, though energy conser- 
vation has been discarded by including all eigenstates of the Hamiltonian in 
the summation, it has been effectively reinstated by the upper limit of the 
Q-integration. 

With the Dirac sea of electrons, we are not so lucky. For each negative 
energy state of momentum pj, there is a corresponding positive energy state 
of momentum pj + q, and the energy difference between these states can be as 
large as 2\Ei\ or larger. Since summation over the entire Dirac sea includes 
the limit — > oo, transitions involving the production of arbitrarily high 
energy electron-positron pairs make significant contributions to the sum. Al- 
though the corresponding matrix elements decrease with increasing \Ei\, they 
do not decrease fast enough to keep the sum over electrons from diverging. 
An explicit calculation, readily done for the free-electron model, shows that 
the resulting sum is infinite. 

Though these arguments have specifically considered a free-electron model, 
the basic ideas hold for bound electrons, as well. The transition probability 
will in this case be broadened by the binding energy, but pairs of arbitrarily 
high energy will still be produced for any given recoil energy Q. Therefore, 
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the conclusions of the previous paragraphs should be valid generally. The 
same conclusions also apply to the TRK sum, since it may be obtained di- 
rectly from the Bethe sum, as shown in Eq. ()14p . In the photoabsorption 
sum of GGT the momentum transfer is fixed by the relationship, q = Eni/c. 
While this differs from that in the free-electron model for charged particle 
scattering, the above arguments do not depend on the precise relationship 
between these quantities. Therefore, the arguments are applicable to all sum 
rules considered in this article. 

I wish to stress that it is completely "reasonable" that the result is infi- 
nite. The use of closure has removed any limits on the energy transfer, and 
the model no longer obeys conservation of energy, at least not with respect 
to the finite energy of any realistic physical system. This point is of some 
interest, since it seems to imply that the presence of infinities is an unavoid- 
able consequence of the use of closure along with inclusion of pair production 
processes in these sums. The sums include contributions from creation of ar- 
bitrarily high energy electron-positron pairs. Furthermore, these are real, as 
opposed to virtual, pairs. It is not, for example, a matter of an electron being 
surrounded by a sea of virtual pairs, leading to polarization of the vacuum 
and altering observable quantities, such as the electron's charge and mass. 
Rather, these are real pairs being produced by collision with a projectile, 
which in this treatment is effectively being considered to have infinite en- 
ergy. Since there is infinite energy available, real pairs of unlimited energy 
can be produced, and the infinite result for the mean energy transfer is in 
this sense quite reasonable. Though one might have hoped that an approach 
utilizing the full power of quantum electrodynamics, including the removal of 
infinities by the methods of renormalization, would allow the effects of pair 
production to be incorporated along with the usual excitation of positive en- 
ergy electrons, these observations seem to indicate that any such efforts will 
prove unsuccessful if they also utilize closure to obtain sum rules. It appears 
that calculations of the pair production contributions to these sums will re- 
quire an approach with explicit limits on the amount of energy exchanged 
between projectile and target. 

6 Conclusion 

In this article, I have discussed relativistic effects on sum rules appearing in 
the description of scattering of charged particles and photons from atoms. 
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Having begun with considerations of how sum rules arise in these contexts, 
I then reviewed the early work on this challenging subject. Following this, I 
discussed my own recent contributions utilizing the Foldy-Wouthuysen and 
Douglas-KroU transformations, and conjectured the latter has the capability 
of producing relatively simple expressions that give highly accurate results 
for these sums for even the heaviest of elements. I also demonstrated the 
equivalence of sums for the transverse interaction in the charged particle 
case and those which describe photoabsorption. Finally, 1 have discussed the 
difficulties that arise in evaluating these sum rules from the perspective of 
Dirac's hole theory, and shown that the understanding of these difficulties is 
much less straightforward then when considered from a single-particle point 
of view. In hole theory, the difficulties were traced to infinities that arise due 
to the presence of the filled Dirac sea. These considerations suggest that it 
will not be possible to obtain physically reasonable sum rules using closure 
if pair production processes are also included in the treatment. 
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